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We present an update of previous work on charged lepton flavor violation 
(LFV) in the seesaw model. The most recent neutrino data fits and post 
WMAP mSUGRA benchmark scenarios are used as input. In this framework 
we compare the sensitivity of rare radiative decays on the right-handed Ma- 
jorana mass scale Mr with the reach in slepton-pair production at a future 
linear collider. 



1 Introduction 

The evidence for flavor violating neutrino oscillations provides strong moti- 
vation to search for LFV also in the charged lepton sector. While charged 
lepton-flavor violating processes are suppressed in the Standard Model with 
right-handed neutrinos [1] due to the light neutrino masses, in supersymmetric 
models new sources of lepton flavor violation exist. For example, the massive 
neutrinos affect the renormalization group equations of the slepton masses 
and the trilincar couplings, and give rise to non-diagonal matrix elements in- 
ducing LFV [2] . This finding has initiated a prospering research activity both 
on LFV in rare decays [3, 4, 5] as well as on lepton-flavor violating processes 
at future colliders [6] . Recently, wc performed a comparative study of the sen- 
sitivity of rare radiative decays on the right-handed Majorana mass scale Mr 
[7] and the reach in slepton-pair production at a future linear collider [8] in the 
context of the SUSY seesaw model. We used mSUGRA benchmark scenarios 
[12] designed for linear collider studies and paid particular attention to the 
uncertainties in the neutrino input parameters. Here we present an update of 
these works, implementing the most recent neutrino data fits and refined post 
WMAP mSUGRA benchmark scenarios. 
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2 Neutrino parameters 



In the last decade a rather unique picture of neutrino mixing has emerged. As 
can be seen in Fig. 1, large to maximal mixing has been established for solar 
and atmospheric neutrinos, while the third angle is strongly constrained by 
reactor measurements. Recently, this picture of neutrino mixing has been re- 
fined further. The results of the KamLAND reactor experiment confirmed the 
disappearance of solar electron neutrinos, while the SNO experiment allowed 
for the first time to study both solar neutrino appearance and disappearance 
via the separate measurement of charged current, elastic scattering and neu- 
tral current interactions. Moreover, the first data of the K2K long baseline 
accelerator experiment indicate a confirmation of atmospheric neutrino os- 
cillations. At the time when a linear collider will be in operation, even more 
precise measurements of the neutrino parameters will be available than today. 
In order to simulate the expected improvement, we take the central values of 
the mass squared differences Amfj = |m| — m|| and mixing angles 6ij from the 
most recent global fit to existing data [10] with errors that indicate the antic- 
ipated 2a intervals of running and proposed experiments as further explained 
in [7]: 

tan 2 23 = l.lO±^, tan 2 13 = 0.006±g;ggJ, tan 2 12 = 0.43±g;& (1) 
Am\ 2 = 6.9±8:|« ■ KT 5 eV 2 , Amj 3 = 2.6±i;| • KT 3 cV 2 . (2) 

Furthermore, for the lightest neutrino mass we assume mi < 0.03 eV, where 
mi ~ eV corresponds to the case of a hierarchical spectrum, while mi ~ 
0.03 eV approaches the degenerate case. 

3 mSUGRA benchmark scenarios 

For our numerical investigations we focus on the mSUGRA benchmark scenar- 
ios proposed in [11] for linear collider studies. These models are consistent with 
direct SUSY searches, Higgs searches, b — ► 57, and astrophysical constraints. 
Recently, the precision measurement of the cosmic microwave background by 
the WMAP experiment combined with previously available data provided a 
refined estimate of the density of cold dark matter in the universe. Assuming 
that most of the dark matter is composed of neutralino LSPs, the smaller 
neutralino relic density led to a shift of previous benchmark points [12] to- 
wards lower values of the universal scalar mass m . In table 1 we specify the 
mSUGRA parameters of the benchmark scenarios B', C, G', and P. These 
are the only models of [11] with left-handed slcptons which are light enough 
to be pair-produced at e + e~ colliders with cm. energies \fs = 500-^800 GeV. 
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Fig. 1. Summary of evidences for neutrino oscillations [9]. 



Scenario 


mi/ 2 /GeV 


mo/GeV 


tan/3 


m 6 /GeV 


m^o /GeV 


B' 


250 


60 


10 


192 


98 


C 


400 


85 


10 


291 


163 


G' 


375 


115 


20 


291 


153 


r 


350 


175 


35 


310 


143 



Table 1. Parameters of selected mSUGRA post WMAP benchmark scenarios [11]. 
The sign of /i is chosen to be positive and Ao is set to zero. Given are also the masses 
of the heaviest charged slepton and the lightest neutralino, that is the LSP. 



4 SUSY seesaw mechanism and slepton mass matrix 

If three right-handed neutrino singlet fields vr are added to the MSSM particle 
content, one gets the following additional terms in the superpotential [4]: 

Wu = -^Mv c R + vfY v L ■ H 2 . (3) 
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Here, Y v is the matrix of neutrino Yukawa couplings, M is the right-handed 
neutrino Majorana mass matrix, and L and H 2 denote the left-handed lepton 
and hyperchargc +1/2 Higgs doublets, respectively. At energies much below 
the mass scale of the right-handed neutrinos, W v leads to the mass matrix 

M v = mF D M- x m D = Y^M^Y^v sin P) 2 , (4) 

for the left-handed neutrinos. Thus, light neutrino masses are naturally ob- 
tained if the typical scale of the Majorana mass matrix M is much larger than 
the scale of the Dirac mass matrix bid = Y^H®), where {H®) = vsin/3 is the 

(H°) 

appropriate Higgs v.e.v. with v = 174 GcV and tan/3 = jj^j- Diagonalization 
of M v by the unitary MNS matrix 

U = diag(e* , , 1)V(0 W , 0i 3 , #23, S), (5) 

where ifii^ and S are the Majorana and Dirac phases, respectively, and 9ij are 
the mixing angles, leads to the light neutrino mass eigenvalues m^: 

U T M U U = diag(mi,m2,m3). (6) 

The other neutrino mass eigenstates are too heavy to be observed directly. 
However, they give rise to virtual corrections to the slepton mass matrices that 
can be responsible for observable lepton-flavor violating effects. In particular, 
the 6x6 mass matrix of the charged sleptons is given by 

/ TO? (TO? ) t \ 



m~ m~ 



with 



{mfjab = {m 2 L ) ab + 5ab(m 2 i a +m|cos2/3 ^-i + sin 2 9 w ^j ^ (8) 

(m~ R ) a b = (m R ) ab + 5 ab (mf a - m% cos 2/3 sin 2 6 W ) (9) 
(mf L J ab = A ab v cos P - 5 ab mi a n tan (3. (10) 



When to? is renormalized from the GUT scale Mx to the electroweak scale 
one obtains, in mSUGRA, 

m| = TOqI + (5m 2 L ) MSSM + 5m 2 L (11) 
m R = m o 1 + {Sm 2 R ) MSS M + Sm 2 R (12) 
A = A Yi + 6A MSSM + SA, (13) 

where mo is the common soft SUSY-breaking scalar mass and A the common 
trilincar coupling. The terms {8m 2 L fl ) M ssM and <L4 MSSM are well-known flavor- 
diagonal corrections. In addition, the evolution generates off-diagonal terms 
in 8m 2 L R and 5 A 2 which in leading-log approximation are given by [5] 
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Fig. 2. Diagrams for 



in the MSSM 



8m 2 L 



- — (3ml + Al)(YjLY u ) 



16tt 2 



with 



(14) 
(15) 
(16) 

(17) 



and Mj, i = 1,2,3 being the eigenvalues of the Majorana mass matrix M. 
In the above, we have chosen a basis in which the charged lepton Yukawa 
couplings and M are diagonal. 

The product of the neutrino Yukawa couplings YJ LY V entering these cor- 
rections can be determined as follows [4]. By inverting (6), one obtains 



Y v = 



v sin P 



R diag (y/mi,y / m2 1 \/m^) W, (18) 



where R is an unknown complex orthogonal matrix. For real R and degenerate 
right-handed Majorana masses, R as well as <f>\ and cf>2 drop out from the 
product YJLY V . Using then the neutrino data sketched in section 2 as input 
the result is evolved to the unification scale Mx- In what follows we refer 
to this illustrative case which suffices for the present discussion. For small 
neutrino masses, mf -C Amf^ the above procedure yields 



WLY v ) a 



M f . 



b - 2 • 2 a I vH^^i 
tr sin p 



1,2 



Aml 3 U a3 U^ 



Upon diagonalization, the flavor off-diagonal corrections in (11)-(13) generate 
flavor- violating couplings of the slepton mass eigenstates. 
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Fig. 3. Diagrams for e+e -> /+Z a -> x a Xf> 



5 Charged lepton flavor violation 



At low energies, the flavor off-diagonal correction (14) induces the radiative 
decays h — > Ijj. From the photon penguin diagrams shown in Fig. 2 with 
charginos/sneutrinos or neutralinos/charged sleptons in the loop, one derives 
the decay rates [4, 5] 



r(k -> Ijj) cx a 3 ml - — r^- tan 2 0, 



(20) 



where m characterizes the sparticle masses in the loop. Because of the domi- 
nance of the penguin contributions, the process [i — > 3e, and also fj,-e conver- 
sion in nuclei is directly related to /i — > ej, e.g., 



Br(fi -> 3e) 
Br([i — > e7) 



^8 m|_ll 
8tt3 I m 2 4 



(21) 



in 



At high energies, a feasible test of LFV is provided by the process e + e~ — ► 
7 — > lj^l^XaX% From the Feynman graphs displayed in Fig. 3, one can 
see that LFV can occur in production and decay vertices. For sufficiently 
narrow slepton widths Fr and degenerate masses, the cross-section can be 
approximated by 



<r(e 



1+ ll)Br{l+ - i+ xo)Br(/^ - ir *„). (22) 



where er(e + e _ — > it IJ) can be replaced by the flavor-diagonal cross-section 
for slepton pair production. The flavor change is described by the factor in 
front of the r.h.s. of (22) resulting in a single mass insertion. 

In the following numerical study we have not assumed slepton degeneracy 
and have summed the amplitudes for the complete 2^4 processes coherently 
over the intermediate slepton mass eigenstates. 
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Fig. 4. Branching ratio of r — »■ /X7 (upper band) and /i — > e7 (lower band) in the 
mSUGRA scenario B'. 
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Fig. 5. Cross-sections for e + e — > /^ + e + 2x? (upper band) and e + e — * r + /i + 
2x? (lower band) at s/a = 500 GeV for the mSUGRA scenario B'. 

6 Numerical Results 

Fig. 4 illustrates the dependence of Br(fi — ► e7) and i?r(r — * /J7) on the 
Majorana mass Mr. The spread of the predictions reflects the uncertainties in 
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Fig. 6. Correlation between Br(fi — > e~f) and <r(e + e — > /i + e + 2x?) at -^s = 800 
GeV for the mSUGRA scenarios (from left to right) C, G\ B' and F. 



the neutrino data. The update of neutrino and SUSY parameters leads to only 
a slight decrease of Br(fi — > ej) as compared to the previous results published 
in [7] . One sees that a branching ratio in the range between the present bound 
and the sensitivity limit of the new PSI experiment, that is 10 -11 <; Br(n — > 
ej) ^ 10 -13 , would point at a value of Mr between 5 • 10 12 GeV and 5 • 
10 14 GeV. On the other hand, Br(r — ► /ij) is more strongly affected by 
the smaller value of mo in the new benchmark models [11], resulting in a 
reduction by a factor of about 5 as compared to the results in [7]. Even if 
Br{r ni) = 1CT 8 is reached, the goal of SUPERKEKB and LHC [13], one 
would only probe Mr <; 10 15 GeV [7]. Nevertheless it is interesting to note 
that r — > /ij is much less affected by the neutrino uncertainties than /i — > e^f. 

Analogously, Fig. 5 shows the Mr dependence of the cross-sections for 
e + e~ — > /i + e~ + 2x1 and e + e~ — > r + /i~ + 2\i- In this case, the [ie channel is 
enhanced by both the larger Am\ 2 and the smaller mo in the new parameter 
set, resulting in a cross-section one order of magnitude larger than what was 
found in [8]. For the r/j, final state, on the other hand, the effect of the smaller 
Am\^ is compensated by the enhancement due to the smaller value of mo, so 
that the net effect is negligible. As can be seen, for a sufficiently large Majorana 
mass Mr the LFV cross-sections can reach several fb. The re channel is 
strongly suppressed by the small mixing angle #13 , and therefore more difficult 
to observe. 

The Standard Model background mainly comes from VF-pair production, 
W production via ^-channel photon exchange, and r-pair production. A 10 
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Fig. 7. Correlation between Br(r — > fjrj) and <r(e + e — ► r + p + 2xS) at yfs = 800 
GeV for the mSUGRA scenarios (from left to right) C\ B', G' and I'. 



degree beam pipe cut and cuts on the lepton energy and missing energy reduce 
the SM background cross-sections to less than 30 fb for /ie final states and 
less than 10 fb for r/x final states. If one requires a signal to background 
ratio Sy V S + B = 3, and assumes an integrated luminosity of 1000 fb^ 1 , a 
signal cross-section of 0.1 fb could only afford a background of about 1 fb. 
Whether or not such a low background can be achieved by applying selectron 
selection cuts, for example, on the acoplanarity, lepton polar angle and missing 
transverse momentum has to be studied in a dedicated simulation. For lepton 
flavor conserving processes one has found that the SM background to slepton 
pair production can be reduced to about 2-3 fb at y/s = 500 GeV [14]. 

The MSSM background is dominated by chargino/slepton production with 
a total cross-section of 0.2-5 fb and 2-7 fb for /ie and r/x final states, respec- 
tively, depending on the SUSY scenario and the collider energy. Here, only the 
direct processes are accounted for. However, the MSSM background in the re 
channel can also contribute to the /xe channel via the decay r — > [iv^v T . If T\ 
and are very light, like in scenarios B' and I', this background can be as 
large as 20 fb. On the other hand, such events typically contain two neutrinos 
in addition to the two LSPs which are also present in the signal events. Thus, 
after r decay one has altogether six invisible particles instead of two, which 
may allow to eliminate also this particularly dangerous MSSM background by 
cutting on various distributions. But also here this needs to be studied in a 
careful simulation. 
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Particularly interesting and useful are the correlations between LFV in 
radiative decays and slepton pair production. Such a correlation is illustrated 
in Fig. 6 for e+e~ — > /i+e~ +2x? and Br(fi — > e7). One sees that the neutrino 
uncertainties drop out, while the sensitivity to the mSUGRA parameters re- 
mains. Furthermore, while models C, G' and I' are barely affected by the 
change in the new parameter set as compared to the set used in [8] , in model 
B' (j(e + e~ — > /i + e~ + 2x") for a given Br{\i — ► e7) is by a factor 10 larger than 
in the previous benchmark point B. An observation of \i — > ej with a branch- 
ing ratio smaller than 1CP 11 would thus be compatible with a cross-section as 
large as 1 fb for e+e~ — ► J2b a K^a ~~ * M + e~ +2x? , at least in model C, G' and 
B'. On the other hand, no signal at the future PSI sensitivity of 1CP 13 would 
constrain this channel to less than 0.1 fb. The correlation of Br(r — > /ry) and 
a(e + e~ — > T + /i~ + 2%?) is shown in Fig. 7. _Br(r — > /ry) < 3 • 10~ 7 does 
not rule out cross-sections in the t/i channel of 1 fb and larger. However, one 
has to keep in mind that these correlations depend very much on the SUSY 
scenario. 



7 Conclusions 

SUSY seesaw models leading to the observed neutrino masses and mixings can 
be tested by lepton- flavor violating processes involving charged lcptons. We 
have presented an updated analysis of the prospects for radiative rare decays 
h ~ * lj"f and slepton pair production and decay e + e~ — ► l^l^ — > +$■ 
Assuming the most recent global fits to neutrino oscillation experiments [10] 
we have illustrated the impact of the uncertainties in the neutrino parameters. 
Furthermore, using post-WMAP mSUGRA scenarios [11] we have investigated 
the dependence of LFV signals on the mSUGRA parameters. For scenario 
B' our results can be summarized as follows. A measurement of Br(fi — > 
e'y) w 10~ 13 would probe Mr in the range 5 • 10 12 -j- 5 • 10 13 GeV, while 
a measurement of Br(r — > ^7) w 10~ 8 would allow to determine Mr ~ 
10 15 GeV within a factor of 2. Furthermore, Br{\i — ► ej) = 10~ 13 -e- 10 -11 
implies a(e+e- -» ^ + e~ +2x°) = 0.02^-2 fb at ^/s = 800 GeV, while Br(r -» 
/try) = 10~ 8 -7- 3 • 10~ 7 predicts cr(e + e- -> + 2%?) = 1 10 fb, again at 

\/s = 800 GeV. Hence, linear collider searches are nicely complementary to 
searches for rare decays at low energies and at the LHC. 
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